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Previous numerical studies have shown that the "ultimate regime of thermal convection" 
can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers 
are eliminated by replacing the walls with periodic boundary conditions (homogeneous 
Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^^^ and tur- 
bulence intensity as Re ~ Ra ' , where the Rayleigh number Ra indicates the strength 
of the driving force. However, experiments never operate in unbounded domains and it 
is important to understand how confinement might alter the approach to this ultimate 
regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally con- 
fined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the 
ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic 
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boundary layers, we still find Nu ~ Re ~ Ra ' . The system supports exact solutions 
composed of modes of exponentially growing vertical velocity and temperature fields, 
with Ra as the critical parameter determining the properties of these modes. Counterin- 
tuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations 
are susceptible to these solutions which can dominate the dynamics and lead to very high 
and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the 
system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. 
velocity and temperature fields. We also test that physical results become independent of 
the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio 
is increased. 



1. Introduction 

There has been longstanding scientific interest in the heat transfer in Rayleigh-Benard 
convection (RBC) - when a fluid conflned between two plates undergoes thermal con- 
vection due to a temperature difference between the cold top and hotter bottom plate 
(Kadanoff (2001); Siggia (1994); Ahlers et al. (2009); Lohse & Xia (2010)). The dynamics 
of the system depend on the strength of the driving temperature difference (given by the 
Rayleigh number, Ra) and the ratio of the kinematic viscosity to the thermal diffusivity 
(given by the Prandtl number, Pr). Recently, studies of unconfined thermal convection 
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- starting with the numerical work by Lohse & Toschi (2003) - were prompted by the 
relevance to natm^al convection phenomena, as occurs in the Earth's atmosphere (Celani 
et al. (2007)) or in the core of stars (Garaud et al. (2010)), and as a way to test theo- 
retical predictions for the scaling of the heat transfer (given by the Nusselt number, Nu) 
and turbulence intensity (given by the Reynolds number. Re). 

In the limit of high Ra it was proposed that the dynamics would reach an ultimate 
regime of thermal convection, which is dominated by the bulk flow rather than the vis- 
cous and thermal boundary layers (Kraichnan (1962)). Within the Grossmann-Lohse 
theory of thermal convection, this regime where the dynamics does not depend on ther- 
mal plumes is also predicted (Grossmann & Lohse (2000, 2001, 2002)). Under these 
conditions, the Reynolds and Nusselt numbers are predicted to scale like Re '-^ VRa and 
Nu ~' VRa, with different Pr scaling depending on the theory (Kraichnan (1962); Spiegel 
(1971); Grossmann & Lohse (2000, 2001, 2002)). The question of the existence of this 
regime was addressed in simulations of homogeneous Rayleigh-Benard convection in a 
tri-periodic cell (i.e., without confinement) where evidence was found supporting such 
a regime when all of the boundaries (and thus boundary layers) were absent (Lohse & 
Toschi (2003); Calzavarini et al. (2005, 2006)) and confirming (Calzavarini et al. (2005)) 
the Prandtl number dependence as suggested by Grossmann & Lohse (2000, 2001, 2002). 
Exponentially growing solutions at low Ra were discovered analytically by Calzavarini 
et al. (2006) and observed in the simulations, appearing as "elevator modes" composed 
of strong upwards and downwards jets, and coinciding with an increasing heat transfer 
until the modes became unstable and disintegrated. 

Later, experiments performed in long rectangular channels by Gibert et al. (2006, 
2009); Tisserand et al. (2010)) and cylindrical pipes by Cholemari & Arakeri (2009) also 
found scaling laws consistent with those expected for the ultimate regime, despite the 
presence of the side walls which cause an additional anisotropy and drag in the flow. 
They also observed a mean global flow composed of hot upwards- and cool downwards- 
moving columns. Whereas in standard RBC strong thermal gradients occur at the top and 
bottom plates, leaving the temperature field within the bulk flow nearly uniform, in the 
experiments with long cells, there exists a mean, linear temperature gradient throughout 
the bulk. It is this underlying linear gradient which is used to drive the convection in 
homogeneous RBC, and in the axially-periodic convection cell considered here. 

In this work, we numerically and analytically investigate thermal convection in a long 
vertical cylinder, in both the low Ra and high Ra limits. Though the top and bottom 
boundary layers are absent, the side walls are still present as in the experimental situa- 
tions. The objective is to find out whether also in this confined (and presumably more 
stabilizing) geometry exponentially growing solutions (so called elevator-modes) exist, 
how the Nusselt and Reynolds number depend on Ra, and how the results compare to 
recent experiments in confined flow geometries (Perrier et al. (2002); Gibert et al. (2006, 
2009); Cholemari & Arakeri (2009)). We first derive exponentially growing mode solu- 
tions of the governing Boussinesq equations with laterally confined boundary conditions, 
and compare them with results from simulations. Then we increase Ra and find that in- 
teraction between the modes stabilizes the system. This interaction prevents the growth 
of the exponential modes, and allows us to define the global heat transfer and the global 
Reynolds number of the system. 
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2. The axially-homogeneous Rayleigh-Benard system and a class of 
exact analytic solutions 

2.1. Definition of the system 

The system under study is axially-homogeneous thermal convection in a vertical cylinder, 
i.e., in a domain with periodic boundary conditions in the axial direction, but with 
lateral confinement. Standard Rayleigh-Benard convection occurs in a cylinder (or other 
geometry) with a top and bottom plate kept at fixed temperatures, and this temperature 
difference drives the flow. In the homogeneous or periodic case, with no boundaries except 
the radial wall, we drive the flow using a temperature gradient A' = dT/dz, where T 
is an underlying linear temperature profile in z about which fluctuations Q occur. The 
absolute local temperature can then be written as T{r, 9, z, t) = T{z) + Q{r, 6, z, t). This 
linear temperature (or density gradient in the case of Cholemari & Arakeri (2009)) is 
observed in experiments (Gibert et al. (2006, 2009)). The diameter of the cylinder d is the 
only imposed scale, and therefore, we make the hypothesis that d is the relevant length 
scale for the convection (aside from the additional smaller scales as turbulence develops). 
The periodicity length L, which can be combined with the diameter to form the aspect 
ratio r = d/L, is much less relevant because it does not change the forcing of the system. 
We non-dimensionalize the governing Navier-Stokes equations using x = d for lengths, 
the free-fall speed u ~ ^ agA' d"^ for velocities, and T ~ A' d for temperatures. The 
Prandtl number is Pr = z^/k and the Rayleigh number is 

agA' d^ , , 

Ra=-5^ . 2.1) 

VK 

Here the material parameters are the coefficient of thermal expansion a, the kinematic 
viscosity and the thermal diffusivity k; the gravitational acceleration is denoted as g. 
Notice that the same deflnition of the Rayleigh number was chosen in the experimental 
studies by Gibert et al. (2006, 2009) and Tisserand et al. (2010). The more usual Rayleigh 
number based on the height of a closed Rayleigh-Benard cell, Rql — agAT L^{vk)~^ , 
where AT is the temperature difference between the horizontal plates, is linked to the 
present one through the relation Rol — Ra F"*. 

In non-dimensional form, the Boussinesq equations read 

— ^-VP+.[^V^u + Qz 
m V Ra 

— = ^^V2e + t(; (2.2) 

where u = (u,v,w) is the velocity field with components in the r, 0, and z directions, 
respectively, and P is the dimensionless pressure field. 

The velocity field satisfies the no-slip condition at the lateral walls located (in non- 
dimensional units) at r = Text = 1/2 so that 



--1/2 



0. (2.3) 



For the temperature T we consider either a fixed condition at the lateral wall, which in 
term of the fluctuations Q means 

Q\r=l/2 = 0, (2.4) 

or adiabatic conditions at the wall, corresponding to: 



d&/dr\r=i/2 = 0. 



(2.5) 
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The system defined so far models the physical situation of a vertical pipe connected to 
two infinite reservoirs of fluids at different temperatures. However, in real experiments 
homogeneous convective systems are produced by connecting a vertical conduct to two 
closed and finite chambers filled with fluid at different temperatures (see e.g. Gibert et al. 
(2006) and Cholcmari & Arakeri (2009)). Therefore, in order to mimic such a condition 
in the following it will be important to take into account the extra boundary condition 
of no mean vertical momentum and no mean temperature increase in the system 

^27r rText /'27r 

w dr de^Q; / / Q dr de = {); (2.6) 

IQ Jo Jo Jo 

which shall be satisfled at any point z on the vertical axis of the cylinder. Note however 
that this "no mean value" condition does not prevent from having a positive heat flow, 
and indeed for the convective regime - as opposed to the conductive one - we expect 
J^-'CwQ dr d0)^ >O. 

2.2. Exact solutions to governing equations 

As in the tri-periodic case, we look for exact solutions which are translationally invariant 
in z and have u = v = (Calzavarini et al. (2006)). The coupled equations for the 
vertical velocity w{r, 9, t) and temperature 0(r, 0, t) are then: 




, ^ , ^. (2.7) 

An adiabatic side wall boundary condition (Eq. 2.5) is relevant for experiments where 
tanks or cells are insulated at the sides, and will be implemented throughout the remain- 
der of this paper. Appendix A gives the analytic result for the fixed temperature side 
wall boundary condition (Eq. 2.4). The mixed (i.e. adiabatic) boundary condition case is 
more difficult to solve analytically, but solutions can still be written for the special case 
Pr = 1. This case will be described in detail in the following. 

Trying solutions with an exponential growth time dependence (e^*), Eq. 2.7 simplifies 
and can be uncoupled by introducing two new fields X{r, 9) and Y{r^ 9) 

X{r,9) = w{r,9) - e{r,9) 

Y{r,9) =w{r,9) + e{r,9) (2.8) 

which satisfy 

9''x idx 1 d^x 

VRa A + l)X = —^ + -— + —7^, 
or^ r or r^ o9^ 

The azimuthal dependence can be expanded as a combination of sin(n0) and cos{n9) and 
written in terms of the modes X„ and K„ which satisfy the modified Bessel equation; 



+ - (VR^(A™„ + l)r'+ n') X„ = 0, 



+ - (VR^(A,„„ -l)r' + n') Y„ ^ 0. (2.10) 

The m label for the growth rate is necessary because again there are multiple solutions 
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for a given n, Xmn, which satisfy the boundary conditions. The radial dependencies of 
the solutions to Eq. 2.10 are 

X„(r) = xo In (V(A,n„ + l)Rai/4 r) 

Yn{r)^yoIn (V(A™„-l)Rai/V) , (2.11) 

where /„ is the modified Bessel function of the first kind. The modified Bessel functions 
of the second kind, are also solutions to Eq. 2.9 but are not physically acceptable 
because they diverge at r = 0. 

The vertical velocity and temperature fields then have the form 

Wnir,t) = ao e^—* [y„(r) + X„(r)] 

e„(r,t) = ao e^""* [Yn{r) - X„(r)], (2.12) 

where the prefactor oq is unprescribcd. 

The growth rate of each mode Xmn depends on Ra. After applying the boundary 
conditions at the cylinder wall at r = Text — 1/2 (Eqs. 2.3 and 2.5), the solutions 
correspond to the roots of 

/„ (^ 7(1+1) Ra^/") _ /„_i(|y(ATT)Ra^/")+/„+i(^y(ATT)Ra^/') 
A + 1 ■ /„(iy(A^Rai/4) /„_i(iy(X^Rai/4)+/„^^(iy(X^Rai/4)' 

(2.13) 

Even though there are multiple solutions A = A,im for each mode ri, the maximum growth 
rate will dominate over the others as time progresses. 

Intuitively one might guess that the exponential growth cannot be maintained, and 
would be arrested by friction on the side walls. However, the solutions to the above system 
of equations have not ignored viscous drag on the walls. Even though the velocity, and 
gradients at the walls, are growing exponentially, the frictional force does not overtake 
the driving buoyancy force. 

In the periodic RB cell, the global, time-averaged dimensionless heat transfer is given 

as 

where the brackets indicate a volume (surface would be the same) and time average. 
When w ^ e^""* and 8 ~ g'*'""*, the instantaneous heat transfer Nu(t) ^ e^'^""* also 
grows exponentially in time. 

2.2.1. Axisymmetric, dipolar modes and critical Rayleigh number 

The spatial form of solutions corresponding to the zero mode case (n = 0) have an 
axisymmetric profile (see Eq. (2.11)). This means that an increasing amount of fluid is 
transported through the cylinder and that the mean temperature of the system increases 
in time. Obviously these zero-mode axisymmetric solutions have a non-zero mean verti- 
cal flow and non-zero mean temperature fluctuation, therefore they are in conflict with 
the condition (2.6). Although these solutions are non-physical (they correspond to an 
unbounded infinite system) they will be useful to validate our simulations. 

From the above considerations we deduce that the first unstable physical mode is n = 1, 
which corresponds to a dipolar flow proflle, moving upwards in half of the cylinder and 
downwards in the other half. Such dipolar modes were also observed in experiments 
in a cylindrical column (Cholemari & Arakeri (2009)). The critical Rayleigh number 
corresponding to Aii(Rac) = can be computed from Eq. (2.13), which leads to the 
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Ra 




m — 1 



2 



3 



7.66 • 10^ 



n = 
1 
2 



0.81 
0.52 
0.14 



0.022 



3.06 • 10* 




1 
2 
3 



0.89 
0.73 
0.52 
0.26 



0.45 
0.12 



1.23 • 10^ 




1 
2 
3 
4 
5 
6 



0.94 
0.86 
0.74 
0.60 
0.44 
0.25 
0.040 



0.70 
0.52 
0.31 
0.080 



0.28 
0.0041 



Table 1. Positive growth rates Xmn corresponding to active modes, for Ra — 7.66 • 10^ (top), 
Ra = 3.06 ■ lO'' (middle) and Ra = 1.23 • 10^ (bottom) at Pr = 1. 



value Rac = 1087.397. In Table 1 Yfe list the active modes at different Ra numbers above 
Rac- For the specific case Ra = 7.66 • 10'^ and Pr = 1.0, the n = maximum growth rate 
is Aio = 0.81. The first physical dipolar mode has grovifth rate An = 0.52. The growth 
rates and their spatial dependencies are confirmed in the following simulations. 

3. Numerical Method 

Numerical simulations in a cylindrical cell are performed using the code developed 
by Verzicco & Orlandi (1996) and Verzicco & Camussi (2003), modified to satisfy the 
periodic conditions at z = 0, L. Furthermore, to satisfy the condition (2.6) the vertical 
flow and mean temperature fluctuation are subtracted at each time step. A similar pro- 
cedure was adopted in the earlier tri-periodic cell simulations by Lohse & Toschi (2003); 
Calzavarini et al. (2005, 2006). A finite difference scheme on a staggered grid is used to 
directly solve the incompressible continuity, momentum equation, and energy equations 
under the Boussinesq approximation using a fractional step method (Verzicco & Orlandi 
(1996); Verzicco & Camussi (2003); Kim & Moin (1985)). The grid is non-uniform in the 
radial direction, but is uniform in z due to the periodicity, allowing the pressure solver 
to take advantage of trigonometric expansions and making the routine more efficient. In 
the simulations we must fix the axial periodicity length L, which is here just a numerical 
parameter rather than a physical scale for the system dynamics. We expect, and will 
find, that as the aspect ratio F = d/L is reduced (by increasing the periodicity length 
L), the system properties become independent of F. 

All the simulations have been started with a velocity field at rest and the temperature 
field consisting only of a random perturbation of maximum amplitude O.OlA'd. 

The new code implementation with periodic conditions was first checked by direct 
comparison of the predicted growth rates and radial forms of w and O from the analysis 
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Figure 1. Growth rate (a) of vertical velocity field w for n = axisymmetric mode as function 
of non-dimensional time - simulation with F = 1/2 (solid black line) approaches theoretical 
prediction of Aio = 0.81 (dashed grey line). Spatial forms (b) of w and temperature fluctuation 
field & for Ra = 7.66 • 10 , simulation and theory coincide. 




t 



Figure 2. Growth rate (a) of vertical velocity field w for n = 1 dipolar mode as function 
of non-dimensional time - simulation with F = 1/2 (solid black line) approaches theoretical 
prediction of Aio = 0.52 (dashed grey line). Spatial forms (b) of w and temperature fluctuation 
field B for Ra = 7.66 ■ 10 , simulation and theory coincide. 

in the previous section, both for cases where n = and n = 1 modes dominated, along 
with the fixed temperature side wall boundary conditions (Appendix A). For Ra = 7656, 
Pr = 1, and F = 1/2, the axisymmetric and dipolar modes are isolated at early times 
and are shown in Figs. 1 and 2, comparing well to the theoretical predictions. Here, 
because the velocity and temperature fields are observed to be z-independent, we use the 
maximum vertical velocity field calculated at each time step to represent ap e^™"* from 
Eq. 2.12. 

When only a few of the exponentially growing modes are active at low Ra (as for 
the parameters in Table 1), the velocity and temperature gradients at On = 2iT/n grow 
exponentially and thus the resolution unavoidably becomes insufficient at some point. 
However, this is only an issue for low Ra cells because as more modes are available, their 
interaction prevents such a situation from occurring. We find that though the resolution 
can affect the maximum w and Q (and thus, the heat transfer) obtained during the 
oscillations of the system, the growth rates are unaffected. For example, at Ra = 7656, 
Pr = 1, and F = 1/4 where the An mode strongly dominates the system, simulations 
with a resolution of 129 x 65 x 129 (in 9, r, z) reach a maximum vertical velocity of 
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50 100 150 

t 

Figure 3. Exponential growth and stabilization of vertical velocity field w for Ra = 7.66 • 10'^, 
Pr = 1, and F = f/4 as function of non-dimensional time, for resolutions 65 x 33 x 129 (solid 
line) and 129 x 65 x 129 (dashed line). Predicted growth rate of An = 0.52 is shown as dotted 
line. 

Wmax — 1000 before stabilizing, compared to Wmax = 200 for a resolution of 65 x 33 x 129 
(shown in Fig. 3). Even though the early time dynamics depend on the resolution, the 
growth rate and steady-state that is reached are unaffected. The situation is analogous 
to what Calzavarini et al. (2006) had found in the unconfined RB system. 

The consistency of the code was also verified using the relation between the Nusselt 
number and the volume and time-averaged kinetic dissipation rate, = i^jVup. For the 
periodic cell, = Nu Ra Pr~^ dr'^ which is confirmed within few percents as long 
as the spatial resolution of the simulation is adequate (see Table 2). It should be noted 
that for increasing Ra the previous identity is verified less satisfactorily since limitations 
of the computational resources prevented us from running cases on larger grids. 

4. Results 

At low Ra wc observe individual growing modes which, depending on Ra and F, domi- 
nate the system dynamics. Already shown in Figs. 1-3 is the exponential growth at early 
times for simulations with Ra = 7.66 • 10"^, Pr = 1, F = 1/2 and 1/4. After the initial 
growth at F = 1/2, Fig. 5 shows that the system temporarily stabilizes as horizontal 
velocity field fluctuations increase, which destroys the axially-uniform modes of Eq. 2.7. 
This process continues periodically, with the horizontal field growing and decaying. This 
is visualized in Fig. 4 which shows snapshots of the dimensionless vertical velocity in 
a vertical cross section through the axis of the cell, just before, during, and after the 
stabilization process has occurred. Analogous dynamics were observed in the tri-periodic 
homogeneous RB system of Calzavarini et al. (2006), and the process was found to depend 
sensitively on the numerical resolution and time-step size. Here, we have an additional 
parameter, the periodicity length which can affect this stabilization process. If the peri- 
odicity length is increased so that F — 1/4, instead of observing the long-term periodic 
behavior, the system settles and stabilizes after the horizontal fluctuations become of the 
same order as the vertical fluctuations (Fig. 5). This signals that the origin of the physical 
process which causes the destruction of the exponentially growing modes, is related to 
the ability of the fluctuations to extend along the z-direction. 

Plotting i(iln(Nu)/(it as a function of time, along with the active Am„, clearly shows 
convergence of the system to the fastest growing mode for Ra = 7.66 • lO'^ with F = 1/2 
in Fig. 6. 
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(a) t = 99.90 



(b) t = 99.96 



(c) t = 100.02 



(d) t = 100.2 
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Figure 4. Instantaneous contour plots of the dimensionless vertical velocity w, through a vertical 
cross section of the cell, for Ra = 7.66 • 10'', Pr = 1, and F — 1/2. Times are chosen just before 
stabilization of the growing dipolar mode a.t t = 99.9 (a), during stabilization when horizontal 
velocity fluctuations are growing at t = 99.96 (b) and t = 100.02 (c), and later after decay at 
t = 100.2 (d). Note that the colorbar scales are varied according to the flow so that the structure 
can be seen at all times. A movie of the process is provided in the Supplementary Material. 




40 60 80 100 



Figure 5. Aspect ratio dependence of mode growth and decay at Ra — 7.66 • 10^, Pr — 1. 
Plotted are the r.m.s. vertical (solid line) and horizontal {uh, dashed line) velocity components 
for F = 1/2 (black) and F = 1/4 (grey). 



As Ra is increased, more modes are present and the growth rates become closer to- 
gether. This leads to interaction between the modes and stabilization of the system by 
preventing individual modes from dominating the flow. Fig. 7 shows the time variation of 
the Nusselt number for Ra = 1.23-10'^ for F = 1/2 and F = 1/4. Though the maximums 
jump among the lower growth-rate modes, there is no prolonged exponential growth as 
was seen at the lower Ra. There is again a difference between F = 1/2 and 1/4 - the 
lower aspect ratio simulations do not reach the same maximum growth rates, leading to 
a more stable system. For reference and comparison to Fig. 4 with F = 1/2, Fig. 8 shows 
a snapshot of the vertical velocity field in a F = 1/8 cell for Ra = 7.85 • 10^. 

We have already seen that the aspect ratio of the cell, determined by the chosen numer- 
ical periodicity length, can affect the stabilization of the system by the breakdown of the 
exponentially growing modes (Figs. 5 and 7). Because we are interested in understanding 
the heat transfer and development of turbulence in a long pipe, when going to higher Ra 
one would hope that the results are unaffected by the choice of the periodicity length. 
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Figure 6. Plot of |dln(Nu)/dt, along with the active mode growth rates from Table 1 for 
Ra = 7.66 • 10^ and Pr = 1, for F = 1/2. The axisymmetric mode with Aio — 0.81 was 
suppressed, so the system periodically sets in on the next fastest growing mode An = 0.52. 



1 
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Figure 7. Plot of idln(Nu)/dt for Ra = 1.23 • 10^ and Pr = 1, for F = 1/2 (panel a) and 
F = 1/4 (panle b). The active mode growth rates reached by the system are shown by the 
horizontal lines at A = 0.60, 0.52, 0.44, 0.31, 0.28, 0.25, 0.080, 0.04, 0.0041 from Table 1. 

Table 2 shows simulation results at Ra ~ 0(10'* — 10®) at aspect ratios F — 1/8, 1/4, 
and 1/2. The aspect ratio was changed by doubling the periodicity length, and keeping 
the resolution in z the same. It appears that the global system properties indeed start to 
converge once F = 1/4—1/8 for Ra = 0(10^ — 10^). For higher Ra, there is still aspect 
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Figure 8. Instantaneous contour plot of the dimensionless vertical velocity w, through a 
vertical cross section of the cell, for a typical state at Ra — 7.84 • 10^, Pr — 1, and F = 1/8 

ratio dependence, but numerical limitations prevent us from increasing the periodicity 
length further. Since the F = 1/4 simulations require half the number of grid points with 
respect to longer cells with F = 1/8, higher Ra simulations which already require a finer 
resolution are done with the intermediate aspect ratio F = 1/4. 

The spatial variation of the r.m.s. velocity components plotted in Fig. 9 for various Ra 
shows the horizontal velocity it/i,rms peaks at the axis, and goes to zero at the walls. The 
vertical velocity fluctuations peak closer to the walls, with a local minimum along the 
axis, falling to zero at the walls owing to friction. At lower Ra, the horizontal fluctuations 
do not reach the same size as the vertical fluctuations at the axis. As Ra is increased, 
u-h.ims increases relative to Wrms at the axis, eventually overtaking Wrms- The horizontal 
fluctuations also become more uniform across the cell. The peaks in Wims become sharper 
and move towards the outside. 

Similar behavior was observed in the measured profiles in the buoyancy-driven pipe 
turbulence experiments of Cholemari & Arakeri (2009). In their experiments with Ra « 
10* and Pr = 670, the profiles resembled those of Fig. 9(a) for the lower Ra = 3.06 • 10^ 
and Pr = 1, but the horizontal fluctuations reached only about 60% the magnitude of 
the vertical fluctuations at the axis. The reason for the resemblance to the lower Ra 
result here is almost certainly the high value of Pr. To check, we also calculated that 
the turbulent shear stress (uhw) was absent across the cell. For all cases, (uhw) is close 
to zero compared to the r.m.s. components. As in the experiments, the system remains 
anisotropic, evidenced by the fact that (w^) / (ufi) 7^ 1, especially near the walls where 
the horizontal fluctuations fall to zero much more rapidly than the vertical fluctuations. 

Countcrintuitively, the heat transfer at low Ra — 0(10"*) and smaller, is very large due 
to single exponentially modes dominating the system in this regime. For these systems, 
it is inappropriate to define an average Nu or Re because we have seen that the maxi- 
mum peaks (as in Fig. 5) depend on the grid resolution. As Ra is increased, the system 

1 /2 

stabilizes due to mode interaction and the scaling appears consistent with Nu Ra ' 
for Ra = 0(10^ — 10*), as shown in Fig. 10. The values used to determine the scal- 
ing law come from the aspect ratio 1/4 runs (given in Table 2). Errors are estimated 
by the percentage difference of the two methods to calculate the Nusselt number (one, 
being the volume average in Eq. 2.14, the other being via the global relation with the 
kinetic energy dissipation rate). The highest Ra runs appear under-resolved but never- 
theless the prediction confirms the scaling laws. The mean Reynolds number is calculated 
as Re = {(u^^y ^Y^"^ d/ V (where u is the total dimensional velocity), and also appears 
consistent with the 1/2 power law from Fig. 10. 

5. Conclusions 



We have presented theoretical results and simulations of thermal convection in a lat- 
erally confined cylindrical pipe in the limit of small aspect ratios. At low Ra we find 
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exponentially growing modes consisting of upwards and downwards-flowing columns, 
analogous to those found in homogeneous RBC with no boundaries and simulations 
in tri-periodic cells (Calzavarini et al. (2005, 2006)), along with experimental obser- 
vations (Gibert et al. (2006, 2009); Cholemari & Arakeri (2009)). The breakdown of 
these modes is found to depend on the numerical aspect ratio of the cell (or periodicity 
length), implying that their destruction is related to the ability of disturbances in the 
axial direction to grow. At higher Ra, the scaling of the heat transfer and turbulence are 
consistent with the predicted (Grossmann & Lohse (2000, 2001, 2002); Lohsc & Toschi 
(2003)) Nu ~ VRa, Re ^ -v/lla scaling of the ultimate, bulk-dominated regime of thermal 
convection, just as in the case of the tri-periodic boundary conditions (Lohse & Toschi 
(2003); Calzavarini et al. (2005)). 
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Figure 10. (a) Scaling of Nusselt number Nu and (b) Reynolds number Re with Ra for F — l/A 
data. Dashed lines have slope 1/2 for comparison. The power law fits (solid lines) have slopes 
of 0.50 ± 0.06 for Nu and 0.51 ± 0.05 for Re. 
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Table 2. Parameters and results from simulations. Columns contain Rayleigh number Ra, aspect 
ratio r, grid resolution (number of nodes in 6, r, and z directions), Nusselt number Nu calculated 
by volume average of Eq. 2.14, Reynolds number Re calculated from mean fluctuating velocity, 
ratio of the vertical to horizontal energies in fluctuating flelds, check of global relation of eu(ViMj) 
calculated directly from velocity gradients compared to eu{Nu) = (u'^ Nu Ra PT~^d~'^), and 
flnally the non-dimensional Kolmogorov scale rj calculated from the kinetic energy dissipation 
rate eu(ViUj). 



It is remarkable that the sidewall boundary conditions and the resulting kinematic 
boundary layers do not lead to different scaling behavior with a less steep increase of Nu 
with Ra, as common for the lower Ra regimes, say Ra < 10^^, in the unifying theory of 
Grossmann & Lohse (2000, 2001, 2002). The reason must be that these kinetic boundary 
layers only form at the sidewalls, and not at the top and bottom walls, which do not 
exist in our setup. We therefore neither expect a less steep increase of Nu with Ra 
beyond the onset of turbulence in the lateral kinetic boundary layers, which is expected 
to occur beyond a critical shear Reynolds number i?e* « 0.5^/Re « 420 (Landau & 
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Lifshitz (1987); Grossmann & Lohse (2000, 2002, 2011)). We note that here we are still 
considerably below this transition, as even for our largest Ra we only have Reg ~ 75 <^ 
Re* . It would be worth-while to push the numerical simulations in this laterally confined 
geometry to Ra numbers beyond the onset of the shear instability, in order to confirm 
that also then Nu ^ Re ^ Ra^^^, without any logarithmic corrections, which seem to 
be typical for ultimate RB flow (and analogous ultimate Taylor- Couette flow) in fully 
confined (i.e., with boundaries and boundary layers also towards the top and bottom) 
geometries (Chavanne et al. (1997, 2001); DubruUe (2001); Grossmann & Lohse (2011); 
van Gils et al. (2011)). 



Appendix A. Zero temperature fluctuations at side wall solution 

Though this boundary condition is not currently relevant for the experiments, it has 
a simple solution which was used to check the numerical results. Zero temperature fluc- 
tuations at the side wall means that physically, the absolute temperature follows the 
imposed linear profile at the walls. Eq. 2.7 is exactly satisfied by solutions of the form 

Wmn{r,0,t) ^ wq e-^™"* cos(n6') Jn{kmnr) 

e™„(r, 61, t) = Go e^'""* cos(n6l) J„(/c„„r) (A 1) 

where J„ is the nth Bessel function of the first kind (sin(n0) dependence is also accept- 
able). The growth rates \mn and wave numbers kmn are determined via the boundary 
conditions, Eqs. 2.3 and 2.4. Here with w{r — 1/2) = 6(r = 1/2) ~ 0, km,n is related to 
the pran (thc mth root of J„) by 

(A 2) 



Text 1/2 

The dispersion relation determining the growth rate for each mode is given by 

2\/PrRa A„™ = -fc„„,' (1 + Pr) + kran (1 + Pr)' + 4 Pr (Ra - k^J) , (A 3) 

which has the same form as the tri-periodic case (up to a constant factor due to the 
different choices for non-dimensionalization) , but the allowed values of k^n differ because 
of the geometry. The critical Raylcigh number for this case, i.e. for Ai4(Rac) = 0, is 
Rac = (fcii)4 = 3448.964. 
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computing time on SARA in Amsterdam. We thank Valcntina Lavezzo for careful reading 
of the manuscript. We also thank the EU COST Action MP0806. 
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